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Abstract 



We study the effect of Coulomb correlations on the ultrafast optical dynamics 
of small metal particles. We demonstrate that a surface-induced dynamical 
screening of the electron-electron interactions leads to quasiparticle scattering 
with collective surface excitations. In noble-metal nanoparticles, it results in 
an interband resonant scattering of (i-holes with surface plasmons. We show 
that this size-dependent many-body effect manifests itself in the differential 
absorption dynamics for frequencies close to the surface plasmon resonance. 
In particular, our self-consistent calculations reveal a strong frequency depen- 
dence of the relaxation, in agreement with recent femtosecond pump-probe 
experiments. 

Pacs numbers: 36.40. Gk, 36.40.Vz, 61.46. +w, 78.47.+p 
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The absorption of light by metal nanoparticles is dominated by the surface plasmon (SP) 
resonance. In small particles, the strong three-dimensional confinement changes both the 
static and dynamic optical properties. For example, the scattering of single-particle excita- 
tions at the surface leads to SP damping and to a resonance linewidth inversely proportional 
to the nanoparticle size Another static effect is the size-dependent enhancement of the 
non-linear optical susceptibilities |2|||. At the same time, the effect of confinement on dy- 
namical properties is much less understood. [I],f|]. One can mention, for instance, the role of 
electron-electron interactions in the process of cluster fragmentation, the role surface lattice 
modes in providing additional channels for intra-molecular energy relaxation, and the in- 
fluence of the electron and nuclear motion on the superparamagnetic properties of clusters. 
Ultrafast optical spectroscopy allows one to address the dynamical aspects by probing the 
time evolution of the excited states with a resolution shorter than the dephasing or energy 
relaxation times. 

Recently, electron relaxation in nanoparticles attracted much interest p|-|TU||. Ultrafast 
pump-probe experiments |7]U performed on noble-metal particles reported a time-dependent 
spectral broadening of the SP resonance, originating from the pump-induced heating of the 
electron gas. The inelastic carrier scattering was found to be essential for understanding the 
differential transmission lineshape JFJ. Similar to metal films JTT|-p^] , the dynamics revealed 



an initial non-equilibrium stage during which the e-e scattering leads to thermalization of 
the electron distribution, followed by a relaxation to the lattice. However, certain aspects 
of the optical dynamics in nanoparticles were observed to differ significantly from films. In 
particular, experimental studies of small Cu nanoparticles [0] revealed considerably faster 
relaxation above and below the SP resonance. These observations suggests that, in strongly 
confined systems, collective surface excitations should play an important role in the electron 
dynamics. 

In this paper we address the role of collective surface excitations in the electron relaxation 
in metal nanoparticles. We demonstrate that the e-e interactions inside small particles are 
strongly modified due to a surface- induced dynamical screening. This leads to a possibility 
of quasiparticle scattering accompanied by the emission of SP's. We show, in particular, 
that in noble-metal particles, this novel size-dependent many-body mechanism results in a 
resonant scattering of the ci-hole into the conduction band. The latter effect manifests itself 
in a strong frequency dependence of the optical dynamics close to the SP resonance, and 
can be observed experimentally with ultrafast pump-probe spectroscopy. 

We consider spherical metal particles, embedded in a medium with dielectric constant e m , 
with radii R sufficiently small, so that only dipole surface modes can be optically excited. 
The optical properties of such a colloid are determined by the dielectric function, e m + 
3pe m (e — e m )/(e + 2e m ), where e(u) = e'{uj) + ie"{u) is that of a single metal particle and 
p <C 1 is the volume fraction occupied by the nanoparticles. For noble metals, e(u>) includes 
the (complex) interband dielectric function, e^cu), associated with the d-electrons. The 
absorption coefficient is then given by |J 

a ( w ) = _9 p ^ e 3/2 bn _L^ j (1) 
c e sp [LJ) 

with 

e sp {u) = e d {u) - u) 2 Ju{u) + ij s ) + 2e m . (2) 
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Here 7 S ~ vpjR characterizes the SP damping by intraband e-h excitations due to the 
surface scattering, and uo p is the plasmon frequency of the conduction (s-p) electrons. The 
SP frequency, u sp , is determined by e' sp {oj sp ) = 

Equation ([]]) describes the linear absorption by large metal clusters. When the electron 
gas is in the quasiequilibrium with the lattice, the temperature dependence of a(uo) also 
determines the time evolution of the differential absorption measured by pump-probe spec- 
troscopy. This is the case for time delays between the pump and the probe optical pulses 
longer than ~ 1 ps [HJ-13]. Since for electron temperature T much smaller than the Fermi 



energy Ep, the intraband contribution to e sp (uj) is nearly T- independent, the SP dynamics 
comes mainly from the T-dependence of the interband dielectric function e^(u;). The latter 
is typically approximated by its bulk expression calculated for non-interacting electrons [|TJ . 
In such model, the SP dynamics is essentially size- independent. 

Here we point out that in small particles, there is a size-dependent contribution to e^iuj) 
that arises from an interband <i-hole scattering mechanism, absent in bulk metals. This 
many-body effect originates from the fact that the collective surface excitations strongly 
modify the dynamically screened e-e interactions inside nanoparticles, as we demonstrate 
below. 

The screened Coulomb potential at point r arizing from an electron at point r' inside a 
spherical particle, Z7 w (r, r'), is determined by the equation fll4"| 



U u (t, r') = U(r - r') + J dndr 2 C/(r - n)IL(ri, r 2 )U u (r 2 , r') 



(3) 



where U(r — r') is the unscreened Coulomb potential and n w (r 1 ,r 2 ) is the polarization 
operator. The latter includes contributions from the conduction and d-band, as well as from 
the surrounding medium, = IV^ + nf, + U™. For simplicity, we assume here step density 
profiles, rij(r) = fii9(R — r), i — c,d, where n c {fid) is the average density of the conduction 
(d-) electrons. For the high frequencies of interest, can be expanded in yielding to 
lowest order [15 n^(r, ri) = — (l/muj 2 )'V[n c (r)V5(r — ri)]. The <i-band contribution to the 
rhs of Eq. (RI) can be found from the relation 



J dT 2 ILi(T,T 1 )U u (T U T') = VPj(r,r'), 



(4) 



where P^(r, r' 



" / " — — Xd(w)Vf/ L , J (r, r') is the (i-band polarization vector (Pf, = outside the 
nanoparticle) and Xd = {^d — l)/47r is the interband susceptibility. U™ can be obtained from 
a similar relation (with P™ = inside the nanoparticle). After some algebra, Eq. @ takes 
the form 



■ >^(r,r') = f/(r-r') + ' ed 



An 



e 2 n c 



driVr 



1 



r i 



rV 1 9(R-r 1 )U u (T 1 ,I f ) 



(5) 



with iJi = 47re 2 n c /m. We solved Eq. @ by expanding U u in spherical harmonics, Ylm{y) 



The solution reads [JTH] 



r ) e 
— + 



E 



47T 
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Y LM (r)Y* M (r' 



(6) 
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where e(u) = €d{uj) — oj 2 /uj 2 . Here cl(c<j) = e' L (uj) +ie'[(u) is an effective dielectric function, 
whose zeros, e' L {uj) = j^-[e' d (cu) — uj 2 /uj 2 } + j^i^m = determine the frequencies of the 
multipole collective surface excitations |J. Note that e" L {ui) contains a correction due to the 
damping of collective excitations by electron-hole pairs, which can be obtained by adding 
Imfl^ to the rhs of Eq. (||) . 

Equation (|6]) represents a generalization of the plasmon pole approximation for spherical 
particles. The two terms in the rhs describe two distinct contributions. First is the usual 
bulk-like screening of the Coulomb potential. The second, size- dependent contribution de- 
scibes a new effective e-e interaction induced by the surface: the potential of an electron 
inside the nanoparticle excites collective surface modes, which in turn act as image charges 
that interact with the second electron. For frequencies close to that of the SP, it is sufficient 
to retain only the resonant dipole term, L — 1. The corresponding surface-induced e-e 
interaction potential is then 

where e sp (oj) = 3ei(w) is the same as in Eq. ([|). Thus, for sufficiently small R, the dynam- 
ically screened Coulomb potential in the nanoparticle is dominated by the SP pole. This 
leads to the possibility of quasiparticle scattering accompanied by the emission of a SP. 

In particular, consider the SP-mediated scattering of the <i-hole into the conduction band. 
This process is described by the Matsubara self-energy due to the potential (^) |14[ 



^M = -3^£^G^ +i ,). (8) 

Here d aa i = (c,a\r\d,a') = (c,a\p\d,a')/im(E^ — E^,) (we set Ti = 1), where \c, a') (\d, a')) 
and E c a (E^) are the conduction (d-) band eigenstates and eigenenergies, respectively, and 
G c a is the Green function of noninteracting conduction electrons. For a quasicontinuous 
spectrum, the matrix element can be approximated by a bulk-like expression, (c, a|p|rf, ct') = 
8aa'(c\p\d) = 5 aa 'fi, the corrections due to the surface scattering being ~ Rr 1 . After the 
frequency summation, ImS^, takes the form 

Tm^m- 9eV N{E? a -E) + f{E%) 

lm ^ {E) ~ -™W) 2 ^ 3 e^Eg-E) ' (9) 

where E^f = Eg—E^, and f(E) and N(E) are the Fermi and Bose distributions, respectively. 

We see from Eq. (|) that the scattering rate of a <i-hole with energy E%, ■jh(E^) = 
Im£^(-E^), depends strongly on the nanoparticle size. Importantly, it exhibits a peak as 
the interband separation E c ^ approaches the SP frequency u sp . The reason for such a sharp 
energy dependence is that the surface potential (0) induces only vertical (dipole) transitions, 
thus restricting the number of allowed final states. 

Let us turn to the absorption spectrum. Consider an optically excited e-h pair with 
excitation energy u close to the onset of interband transitions, A. As discussed above, the 
d-ho\e can subsequently scatter into the conduction band by emitting a SP. Eq. (^) indicates 
that for u ~ u> sp ~ A, this process will be resonantly enhanced. In particular, this is the case 
in Cu and Au, where u) sp and A are close to each other J7H§. At the same time, the electron 
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scatters in the conduction band via the usual two-quasiparticle process. For u ~ A, the 



electron energy is close to E F , and the scattering rate is j e ~ 10~ 2 eV [ 17| . Using the bulk 
value of n, 2fi 2 /m ~ 1 eV near the L-point [|l!|, we find that 7/1 exceeds 7 e for nanoparticles 
smaller than ~ 5 nm. 

We now calculate the effect of the SP-assisted interband scattering on the absorption 
coefficient a. For non-interacting electrons, the interband susceptibility, Xd{i^) = Xd{i^) + 
Xd{— iu), has the standard form 

= - E zMttk 5: Gt(iu')GUiu> + iuj). (10) 

Since the (i-band is occupied for all energies, it is sufficient to include only the effect of <i-hole 
scattering into x d {iuj) 0. Substituting G d (iuj') = [iu' - E d a + E F - _1 , with E d (iuj) 



given by Eq. (||), we obtain after the frequency summation 

_ eV r dE c g(E c ) f{E c )-l 

Xd[UJ> m 2 J {E cd ) 2 uj - E cd + i lh (u, E c ) ' 1 ' 

where g(E c ) is the density of states in the conduction band, and we assumed a dispersionless 
rf-band with energy E d . Here 7^,(0;, E c ) = lmS d (E c — u) is the scattering rate of a <i-hole 
with energy E c — u>, 

^) = ~lS^ m)lm ^y (12) 

The comparison of Eqs. (|T2|) and (|l|) shows that the scattering rate depends on the frequency 
in the same way as the absorption coefficient: both 7/1(0;, E c ) and a{uj) exhibit a peak at 
uj = uo sp . This implies that the c?-hole experiences resonant scattering into the conduction 
band as the optical field frequency approaches the SP frequency. We thus arrive at a self- 
consistent problem defined by Eqs. ([[]), (0), flTT|), and ([12]). 

The main role of the SP-assisted <i-hole scattering is to change the absorption lineshape 
in the vicinity of the SP resonance. The SP width is determined by e'^(uj) = 47TXd(w), which 
now, according to Eqs. (|T2|) and ([TT|) , acquires a sharp frequency dependence for to ~ u sp . 
It is also important that the effect of 7^ on e'^(uj) increases with temperature. Indeed, the 
Fermi function in the rhs of Eq. ( |12"D implies that 7^ is small unless E c — Ep < k^T ' . Since 
the main contribution to Imx^ es (u;) comes from energies E c — Ep ~ uj — A, the scattering 
becomes efficient for temperatures /c#T > u sp — A. The combined effect of the above uj and 
T dependence manifests itself strongly in the differential absorption dynamics, as illustrated 
below. 

In Fig. 1 we show the results of our self-consistent numerical calculations of the differential 
transmission dynamics. We use the parameters of the experiment , performed on ~ 5 nm 
Cunanoparticles, with the SP frequency, uj sp ~ 2.22 eV, slightly above the onset of the 
interband transitions, A ~ 2.18 eV. Fig. 1(a) shows the spectra at various time delays, 
calculated as the difference between a{oj) at T = 300 K and T(t), the pump-induced hot 
electron temperature. The latter was obtained from the two-temperature model with bulk 
Cu parameters p(J and initial temperature T = 800 K. The d-ho\e scattering leads to 



a steeper cu-dependence of the differential transmission for uj > oj sp ] its effect, however, 



5 



is best seen in the time evolution as described below. We also included the effect of the 
intraband e-e scattering; for electron energy close to Ep, this can be achieved by adding the 
e-e scattering rate 7 e (£ c ) oc [1 - f(E c )][(E c - E F ) 2 + (nk B T) 2 } to -y h in Eq. ([XX]) . The 
difference in ^ e (E c ) for E c smaller and larger than Ep leads to an asymmetric differential 
transmission lineshape, in agreement with experiment ]7| . Note that in small nanoparticles, 
7 e is expected to exceed its bulk value which could account for the more symmetric 
lineshape observed in 40 nm nanoparticles M. 

In Fig. 1(b) we show the calculated time evolution of the differential transmission. The 
relaxation is slowest at uj = uj sp , which can be attributed to the robustness of the SP 
mode; the corresponding relaxation time is 3.5 ps. More importantly, for short time delays, 
the relaxation is significantly faster below and above uj sp , with a characteristic time of 1.4 
ps. These results are consistent with the experimental data shown here for comparison 
in Fig. 1(c) (the experiment also revealed an initial non-equilibrium regime not addressed 
here). In contrast, in the absence of <i-hole scattering [7/, = in Eq. (|TT])], the frequency 
dependence would be smooth above the resonance [see Fig. 1(d)]. Such a smooth dependence 
originates from the lineshape of the absorption peak: since uj sp ~ A, the absorption is 
larger for uj > uj sp due to the interband transitions, and therefore the relative change in the 
absorption is smaller. This, however, changes with ci-hole scattering turned on (7^ 7^ 0). To 
understand this difference, one should note that away from uj sp , the relaxation characterizes 
the temperature-induced change in the SP width. The latter acquires an additional frequency 
dependence near the resonance, originating from that of 7^ (u). This effect is enhanced 
above the resonance due to the similar frequency dependence of 7^ and a. Furtheremore, 
the efficiency of this mechanism increases with T, as discussed above. This results in a 
strong frequency dependence of the relaxation (for uj > uj sp ) for shorter time delays, which 
correspond to higher T [see Fig. 1(b)]. 

In conclusion, we have shown that surface-induced electron-electron interactions in small 
metal particles lead to quasiparticle scattering mediated by collective surface excitations. In 
noble-metal particles, this size-dependent mechanism results in a surface-plasmon-assisted 
resonant scattering of ci-holes into the conduction band. The latter effect can be observed 
with ultrafast pump-probe spectroscopy. Our work points out the increasing role of many- 
body correlations with decreasing size, important for understanding the transition from 
boundary-constrained nanoscale materials to molecular clusters. 

This work was supported by NSF CAREER award ECS-9703453, and, in part, by ONR 
Grant N00014-96-1-1042 and by Hitachi Ltd. 



6 



REFERENCES 

[1] U. Kreibig and M. Vollmer, Optical Properties of Metal Clusters (Springer, Berlin, 1995), 

and references therein. 
[2] F. Hache et al, J. Opt. Soc. Am. B 3 1647 (1986). 
[3] L. Yang et al, J. Opt. Soc. Am. B 11 457 (1994). 

[4] See for instance Physics and Chemistry of Finite Systems: From clusters to Crystals, 
eds. P. Jena et al. (NATO Advanced Study Institute Series C, Kluwer Academic, Dor- 
drecht/Boston, 1992); see also W. A. De Heer, Rev. Mod. Phys. 65, 611 (1993). 

[5] T. Tokizaki et al, Appl. Phys. Lett. 65, 941 (1994). 

[6] T. W. Roberti, B. A. Smith, and J. Z. Zhang, J. Chem. Phys. 102, 3860 (1995). 
[7] J.-Y. Bigot et al, Phys. Rev. Lett., 75, 4702 (1995). 

[8] T. S. Ahmadi, S. L. Logunov, and M. A. Elsayed, J. Phys. Chem. 100, 8053 (1996). 

[9] M. Perner et al, Phys. Rev. Lett. 78, 2192 (1997). 
[10] M. Nisoli et al, Phys. Rev. Lett. 78, 3575 (1997). 
[11] W. S. Fann et al, Phys. Rev. B 46, 13 592 (1992). 
[12] C. K. Sun et al, Phys. Rev. B 50, 15 337 (1994). 

[13] R. H. M. Groeneveld, R. Sprik, and A. Lagendijk, Phys. Rev. B 51, 11 433 (1994). 
[14] G. D. Mahan, Many-Particle Physics (Plenum, New York, 1990). 
[15] A. A. Lushnikov and A. J. Simonov, Z. Phys. 270, 17 (1974). 
[16] T. V. Shahbazyan, I. E. Perakis, unpublished. 

[17] S. Ogawa, H. Nagano, and H. Petek, Phys. Rev. B 55, 10 869 (1997). 
[18] H. Ehrenreich and H. R. Philipp, Phys. Rev. 128, 1622 (1962). 
[19] A. Zangwill and P. Soven, Phys. Rev. B 24, 4121 (1981). 
[20] G. L. Easley, Phys. Rev. B 33, 2144 (1986). 

[21] D. Pines and P. Nozieres, The theory of quantum liquids, (W. A. Benjamin, Inc., New 

York, 1966), Vol. I. 
[22] U. Sivan, Y. Imry, and A. G. Aronov,, Europhys. Lett. 28, 115 (1994). 



7 



FIGURES 



FIG. 1. (a) Calculated differential transmission spectra for positive time delays with initial 
hot electron temperature To = 800 K, e m = 2.25, and 7 S = 0.1 eV. (b) Temporal evolution of the 
differential transmission for frequencies equal, above, and below the SP resonance, (c) Experimental 
temporal evolution of the pump-probe signal Q. (d) Calculated differential transmission in the 
absence of interband <i-hole scattering. 
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